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Abstract 

This paper aims to study a family of Leray-a models with periodic bounbary conditions. 
These models are good approximations for the Navier-Stokes equations. We focus our at- 
tention on the critical value of regularization that garantees the global well-posedness 
for these models. We conjecture that 9 = ^ is the critical value to obtain such results. 
. When alpha goes to zero, we prove that the Leray-a solution, with critical regularization, 

gives rise to a suitable solution to the Navier-Stokes equations. We also introduce an 
interpolating deconvolution operator that depends on "0" . Then we extend our results of 
existence, uniqueness and convergence to a family of regularized magnetohydrodynamics 
C*~) \ equations. 
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1 Introduction 

The dynamics of fluids provide various highly challenging theoretical, as well as experimen- 
tal and computational problems for engineers, physicists and mathematicians. It is widely 
believed that all the informations about turbulence are contained in the dynamics of the 
solutions of the Navier-Stokes equations (NSE) for viscous, incompressible, homogenous 
fluids. 

The NSE for a homogenous incompressible fluid are usually written as 
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O 



(l.l) 



du 

— — h u ■ Via — vAu + Vp = /, 
at 

V-u = 0, 
u t=o = u o- 



- 1—1 

x. 

I The unkowns are the velocity vector field u and the scalar pressure p. The viscosity u, 

the initial velocity vector field uq and the external force / are given. 
The dynamics of several conducting incompressible fluids in presence of a magnetic field are 
described by the magnetohydrodynamics equations (MHD). The MHD involve coupling 
Maxwell's equations governing the magnetic field and the NSE governing the fluid motion. 
The system has the following form 

d t u - ^Au + (t* • V)u - (B ■ V)B + Vp + ^V\B\ 2 = 0, 

(1.2) <! d t B-u 2 AB + (u-V)B-(B-V)u = 0, 

V -B = V-u = 0, 

B t =0 = B , U t=0 = Uq. 
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The unkowns are the velocity vector field u, the magnetic vector field B and the scalar 
pressure p. The kinematic viscosity v±, the magnetic diffusivity 1/2, the initial velocity 
vector field uq and the initial magnetic vector field Bq are given. 

The problem of the global existence and uniqueness of the solutions of the three-dimensional 
Navier-Stokes equations are among the most challenging problems of contemporary math- 
ematics. The first attempt was done to Leray [25] who established the existence of global 
weak dissipative solutions to NSE. Such solutions are called "turbulent" and allow singu- 
larities in the velocity fields. The extension to no-slip boundary condition was done by 
Hopf |18j . Short time regularity of solutions has also been known for many years, as they 
have various interesting partial and conditional regularity results |2U[ [3"5"1 [HI ES] • Histori- 
cally, the partial regularity result of Scheffer [33] proving that the potentially singular set 
in time has zero half-dimensional Hausdorff measure, led to that of Caffarelli, Khon and 
Nirenberg [1] . In [3] they worked with a particular class of weak solutions of the Navier- 
Stokes equations called suitable weak solutions. By a suitable weak solution they mean a 
weak solution of the NSE such that for all T £ (0, +oo] and for all non negative fonction 
(j) E C°° compactly supported in space and time, the following inequality is valid: 

2v I j \Vu\ 2 <j) dxdt < j j \u\ 2 ((p t + vA(f>) + (\u\ 2 u + 2pu) • V</> dxdt 
(13) Jo Jt 3 Jo Jt 3 

+2 1 I fu(f> dxdt 
Jo Jt 3 



The inequality (jl.3p is known as local energy inequality. Caffarelli, Khon and Nirenberg 
[3] proved that for the NSE the singular set of a suitable weak solution has parabolic 
Hausdorff dimension at most equal to 1. This implies that if singularities do exist, they 
must be relatively rare. 

Let us mention that it is not known whether the weak solutions satisfy the local energy 
inequality. Whereas it is remarkable that the weak solutions constructed by Leray [2S] , by 
regularizing the nonlinear term, or by adding hyperviscosity are actually suitable [3j,fTQ^fT6] . 
These problems are also open and very important for the MHD equations |11[I31| although 
the MHD equations are not on the Clay Institute's list of prize problems. 
In order to study various aspects of turbulence, and motivated by the difficulties related 
to the complexity of the 3D Navier-Stokes equations, simpler models and simplifications 
of hydrodynamics equations have been proposed over the years for exemple |29t [2T| [6j [T4"l 

[111 El [191 El ESI EH E21 ET] - 

In this paper, we study a Leray-a approximation for the Navier-Stokes equation subject 
to space-periodic boundary conditions. 
The Leray-a equations we are considering are 



(1.4) 



— + u Vit - vAu + Vp = f inM + xT3, 
ot 

u = a 2e (-A) e u + u, in T 3 

V-u = V'« = 0, / u= / 11 = 0, 

Jr 3 Jt 3 
u(t,x + Lej ) = u(t,x), 

u t=o = u o- 



Where {ej, j=l,2 or 3} is the canonical basis of M 3 , L > 0. 

We consider these equations on the three dimensional torus T3 = (R 3 /73) where T3 = 
2-kI? /L , x € T3, and t 6 [0, +oo[. The unkowns are the velocity vector field u and the 
scalar pressure p. The viscosity v, the initial velocity vector field uq and the external force 
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/ with V • / = are given. 

Let be a non-negative parameter. The nonlocal operator (—A) 61 is defined through the 
Fourier transform 

(1.5) (-Apu(fc) = \k\ 29 u{k). 

Fractionnal order Laplace operator has been used in another a model of turbulence in [33J. 
Existence and uniqness of solutions of other modifications of the Navier-Stokes equations 
have been studied by Ladyzhenskaya [21] Lions [29], Malek et al. [32] . 
Motivated by the above work [33], we study here the case = ^ of (|1.4j) . 
We would like to point out that after having finished writing this paper, we were informed 
that there exists an analysis of a general familly of regularized Navier-Stokes and MHD 
models [T7]. In our paper, we note that the Leray-a family considered is a particular case 
of the general study in [T7] where the results do not recover the critical case = \- The 
regularization critical values of the other a models studied in [T7] will be reported in a 
forthcoming paper. 

Therefore, the initial value problem we consider in particular is: 

— + u- Vtt - vAu + Vn = / inl + xT 3 , 
u = a 2 (— A)±u + u, in T3 
V-it = VU = 0, / u= / u = 0, 

u t=0 = u 0' 

or equivalently (II. 4h with 8 = | and we are working with periodic boundary conditions. 
One of the main results of this paper is to establish the global well-posedness of the so- 
lution to equations (jl.4[) in L 2 (Ts) 3 for 6 = |. We conjecture that 6 = \ is the critical 
value to obtain the above result to eqs. (jl.4p . Therefore, smooth solutions of eqs. (|1.4|) 
with > I do not develop finite-time singularities. 

When a = 0, eqs. (jl,4p are reduced to the usual Navier-Stokes equations for incompress- 
ible fluids. 

We recall that H^iT^) 3 is a scale-invariant space for Navier-Stokes equations, i.e. if 
u(x,t) is a solution to the Navier-Stokes equations, then so is U\(x,t) := Xu(Xx, X 2 t) for 
any A > 0. Kato (|40j.|41j) shows the importance for a functionnal space invariance by 
scaling. A lot of works followed which can be summed up in the following way: if the 
initial data are small in an invariant norm with respect to the scaling then the solution is 
smooth for all time. 

Another result of this paper is to establish the global well-posedness of the solution to eqs. 
(jl.4p in H2(Ts) 3 for & = \ and without smallness conditions on the initial data. 
We also discuss the relation between the Leray-a equations with > \ and the NSE. 
The first result of convergence of a a model to the Navier-Stokes equations is proved in 
[13] where the authors show the existence of a subsequence of weak solutions of the a 
model that converge to a Leray-Hopf weak solution of NSE. We improve their result by 
proving an LP convergence property in order to show the convergence of the weak solution 
to a suitable weak solution to the NSE. Moreover, we extend the theory to other models 
where we consider a deconvolution type regularisation to the NSE and an extension to 
the MHD equations is also given. The deconvolution operator was introduced by Layton 
and Lewandowski [23]. We will define the interpolating deconvolution operator in order 
to regularize the NSE in section fTTTI and the MHD equations in section [8j As the Leray— a 
regularization we will show an LP convergence property to the deconvolution operator, 



(1.6) 
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in order to deduce that the deconvolution regularization will give rise to suitable weak 
solutions. 

The critical values of regularization 9 = j holds for the MHD equations and the Navier- 
Stokes equations so it is natural to ask if the singularities in MHD equations are similar to 
the singularities in the Navier-Stokes equations. A partial answer is given in [J2], where 
an analog of the known Caffarelli, Khon and Nirenberg result is established to the MHD 
equations. The extending of the partial regularity result of Scheffer [M] to the MHD equa- 
tion is studied in [2]. 

The paper is organized as follows. In section 2 we recall some preliminary results used 
later in the proof. In section 3, we present the main existence results (Theorem 13.11 and 
Theorem 13.21 below ). Sections 4 and 5 are devoted to the proof of these two results where 
we use the standard Galerkin approximation. In section 6 we show an LP convergence 
result in order to prove that the Leray— a equations give rise to a suitable solution to 
the Navier-Stokes equations. In section 7 we introduce the interpolating deconvolution 
operator and we give applications to the Navier-Stokes equations. 
Finally, we study a deconvolution regularization to the MHD equations. 



2 Functional setting and preliminaries 

In this section, we introduce some preliminary material and notations which are commonly 
used in the mathematical study of fluids, in particular in the study of the Navier-Stokes 
equations (NSE). For a more detailed discussion of these topics, we refer to [8ll9| [T5 | [26 | [39] . 
We denote by L P (T3) 3 and H s (Ts) 3 the usual Lebesgue and Sobolev spaces over T3. For 
the I? norm and the inner product, we write || • ||^2 and (-,-)l 2 - The norm in H 1 ^^) 3 
is written as || • ||#i and its scalar product as (•, We denote by (H 1 ^^) 3 )' the dual 
space of H 1 ^^) 3 . Note that we have the continuous embeddings 

(2.1) ^(Ts) 3 L 2 (T 3 ) 3 (F 1 ^) 3 )'. 

Moreover, by the Rellich-Kondrachov compactness theorem (see, e.g., |12j). these embed- 
dings are compact. 

Since we work with periodic boundary conditions, we can characterize the divergence free 
spaces by using the Fourier series on the 3D torus T3. We expand the velocity in Fourier 
series as: 

u(t, x) = Uk exp {ik ■ x} , 
kei 

where J = |fc such that k = — — , a € Z 3 , a 7^ ol , 

where for the vanishing space average case, we have the condition 

Uq = u = 0. 

Jt 3 

Since the vector u(t, x) is real-valued, we have 

= u* k for every k, 

where u* k denotes the complex conjugate of u^. 

In the Fourier space, the divergence-free condition is 

■Ufc • k = for all k. 



4 



For s € R, the usual Sobolev spaces iP(T3) 3 with zero space average can be represented 



as 



H s = {u = iik ex P {ik " x }i &-fe = Wjl, ||i*|||fs < oo}, 
kei 



where 



Let 



kei 



V s = {u € if s , u k k = 0}, 
we identify the continuous dual space of V s as V~ s with the pairing given by 

{u,v) V s = ^2 \ k \ 2s Uk ■ V- k . 
kei 

We denote V° by H and V 1 by V, where the norm || ■ |||f = || • ||x,a and the scalar product 
as ("> ")h = (' 5 ')l 2 - The norm in V is written as || • \\y = \\ ■ \\ 2 Hl and its scalar product as 

We denote by P CT the Leray-Helmholtz projection operator and define the Stokes operator 
A := -P a A with domain V(A) := H 2 (T 3 ) 3 n V. For u G we have the norm 

equivalence ||Au|| i 2 = ||w||y2. Furthermore, in our case it is known that A = — A due to 
the periodic boundary conditions (see, e.g., [51 [T3]). 

It is natural to define the powers A s of the Stokes operator in the periodic case as 

A s u = ^2 \k\ 2s Uk exp {ik ■ x} . 
kei 

For u £ U(A S / 2 ), we have the norm equivalence || v4 s / 2 tx. ||^2 = ||ti||v s - 

In particular for s = 1 we recover the space V, hence V{A l l 2 ) = V, (see, e.g., [13]). 

Note also that H s+e is compactly embedded in H s (resp. V s+e is compactly embedded 
in V s ) for any e > 0, and we have the following Sobolev embedding Theorem (see [T]). 

Theorem 2.1 The space H 1 ^ 2 is embedded in L 3 (T3) 3 and the space L 3 / 2 (T 3 ) 3 is embed- 
ded in H" 1 / 2 

The following result deals with some interpolation inequalities. 
Lemma 2.1 Let T > 0, 1 < p\ < p < P2 < oo, s\ < s < S2 and n €]0, 1[ such that 

1 77 1 — 1] 

- = 1 and s = ns\ + (1 — ms2- 

P Pi P2 

2 

Let u € p| L Pi ([0,T],H Si ), then u £ LP([0,T], H s ) and 



i=l 



, ^|| ||?7 II ||1— Tj 

\\ u \\Lp([0,T],H s ) S o ll' u llLPi([0,T],H s i)ll' U llLP2([0,r],/f s 2)- 

TTie result holds true when we work with the spaces L Pi ([0,T],V Si ). 

The regularization effect of the nonlocal operator involved in the relation between u and 
u is described by the following lemma. 
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Lemma 2.2 Let 9 G M + , s£l and assume that u G H s . Then u G H s+2e and 

. . .... 1 . 

(2.2) < -^q\\u\\ H s. 



Proof. When u = exp {ifc • x}, then 



fee/ 



(2.3) M = E i + a 2ei fc |2e explifc-cc}. 

fee/ 1 

Formula (|2.3|) easily yields the estimate (|2.2|) . 

Remark 2.1 For 6 = \, u G H s+ z and I lit II „,i < llull h 3 ■ 

Following the notations of the Navier-Stokes equations, we set 

(2.4) B(u,v) := P a (u-Vv) 
for any u, v in V. 

We list several important properties of B which can be found for example in [15j . 



Lemma 2.3 The operator B defined in (|2.4p is a bilinear form which can be extended as 
a continuous map B : V x V — s> V' . 
For u, v, w G V , 

(2.5) (B(u, v),w) v -i = - (B(u, w), v) v -i , 
and therefore 

(2.6) (B(u,v),v) v -i =0. 
The following result holds true. 



Lemma 2.4 The bilinear form B defined in \2.J$ satisfies the following: 

(i) Assume that u G , v and w G V, Then the following inequality holds 

(2.7) \(B(u,v),w)\ < C||ti|| vi |M|v|M|v- 

(ii) B can be extended as a continuous map B : V? x V? — > V~^. In particular, for 
every u G V% , v G V? the bilinear form B satisfies the following inequalities: 

(2.8) \\B(u, v)\\ i < Cllull i ||i>|| a. 
Proof (i) We have 



(2.9) \(B(u,v),w)\ < 



u (g> v : Vw 

3 



By Holder inequality combined with Sobolev embedding theorem we get 

\(B(u,v),w)\ < \\u ® v \\ L 2 \\w\\ H i 
(2.10) < ||'u|| L 3||'u|| i 6||itf||#i 

< CllitH i llwllvllwllv- 

— II \\ v ^ II II V|| || V 
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(ii) We have 

||-B(u,v)|| 1 < c||tt||r3||V'u||r3 

(2.11) ~ , „ 

< C U l \\v\\ 3 . 

Where he have used Holder inequality and the Sobolev injection given in Theorem 12.11 

The following result is given for any u G L°°([0, T];H) n L 2 ([0, T]; V). 

Lemma 2.5 Ifu G L°°([0, T]; J/)nL 2 ([0, T]; V), thenu G L°°([0, T]; V^)nL 2 ([0, T]; V"i) 
andB(u,u) G L 2 ([0, T]; F" 1 ). 

Proof. The regularity on It is obtained from Lemma 
Then we use Lemma 12^41 to obtain 



Now, 



(B(«,ti))|| v -i <C\\u\\ vi \\u\\ v 



\u\\ i G L°°(0,T) and \\u\\ v G L 2 (0, T), 



implies that ||u || i||«||v G L 2 (0, T) (that is L°°([0, T]; V^)nL 2 ([0, T]; V) C L 2 {[0,T}; H) 

)• 

We remove the pressure from further consideration by projecting (|1.4h by P a and searching 
for solutions in the space V s . 

Thus, we obtain the infinite dimensional evolution equation 



(2.12) 



du 

— + vAu + B(u,u) = /, 
u = a 2e A e u + u. 



. u t=o — u o- 

Note that the pressure may be reconstructed from u and u by solving the elliptic equation 

Ap = V ■ («• Vu). 

Later in Sect. 4 and Sect. 5, we will prove that eqs. (|1.4p have unique regular solutions 
(u,p) such that the velocity part u solve the eqs. (|2.12|) . 

We conclude this section with the interpolation lemma of Lions and Magenes [30] and a 
compactness Theorem see [36]. 

Lemma 2.6 Let s > and suppose that 

fill 

ueL 2 ([0,T];V s ) and — G L 2 ([0, T], V~ s ). 

Then 

(i) u G C([0,T];H), with 

( du 
sup \\u{t)\\ H < c I ||w||l2( [0)T ];v-) + I|-HIl2([o,t],v- s ) 



tG[o,r 

and 

/■ m ^ || i|2 2 / ^ u 
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To set more regularity for the unique solution, we also need the following similar result. 
Corollary 2.1 Assume that for some k > and s > 0, 

u G L 2 ([0, T]; V k+S ) and ^ G L 2 ([0,T],V k ~ s ). 

Then 

(i) u G C{[0,T};V k ), and 

The result holds true in the space H s . 



Theorem 2.2 Let X CC H C Y be Banach spaces, where X is reflexive. Suppose that 
u n is a sequence such that it is uniformly bounded in L 2 ([0,T]; X), and — is uni- 
formly bounded in L 2 ([0,T],Y). Then there is a subsequence which converges strongly 
in L 2 ([0,T],H). 

The above Theorem is a special case of the following more general result: 

Theorem 2.3 (Aubin-Lions) Let X CC H C Y be Banach spaces, where X is reflexive. 
Suppose that u n is a sequence such that it is uniformly bounded in L p ([0, T\; X), (I < p < 
du n 

oo), and ^ is uniformly bounded in L q ([0, T], Y), (1 < q < oo). Then u n is relatively 
compact in L p ([0, T], H). 

3 Main existence theorems 

In the following, we will assume that a > 0, T > 0. 

One of the aims of this paper is to establish the global well-posedness of the solution to 
eqs. (|1.4p in L 2 (Ts) 3 for 9 = -r. It must be mentioned that the result holds true for 6 >\. 
We conjecture that 6 = ^ is the critical value to obtain the following result to eqs. (|1.4p . 

Theorem 3.1 For any T > 0, let f G L 2 ([0, T], V" 1 ) and u G H. Assume that 6 = \. 
Then there exists a unique solution (u,p) := (u a ,p a ) to eqs. jl-4\ ) that satisfies u G 

C([0,T];H)nL 2 ([0,T];V) and ^ G L 2 ! ([0,T];V - 1 ) and p G L 2 ([0, T], L 2 (T 3 )). 
Such that u verifies 



du ,_ 

— + vAu + B(u, u) — f, i 



for every <fi G V and almost every t G (0, T). Moreover, this solution depends continuously 
on the initial data Uq in the L 2 norm. 

Therefore, smooth solutions of eqs. (]1.4p with 6 > j do not develop finite-time singulari- 
ties. 

We also prove the global well-posedness of the solution to eqs. (jl.4p in H 2 (T3) 3 for = \ 
and without smallness conditions on the initial data. 
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Theorem 3.2 For any T > 0, let f € L 2 ([0, T], V~?) and u € Vs. A ssume that 
6 = \. Then there exists a unique solution (u,p) := (u a ,p a ) to eqs. fil.4\) that satisfies 

u eC([0,nv'2)nL 2 ([0,T};Vl), ^ e L 2 ([0,T],V-h) andp£L 2 ([0,T],H-2(T 3 )). 
Such that u satisfies 

^- + vAu + B(u,u)- f,(f>) =0 

at I y-2,V2 

for every <p € V? and almost every t € (0, T). Moreover, this solution depends continu- 
ously on the initial data uq in the norm. 

We note that the above result holds true for any s > and without smallness conditions 
on the initial data. 



4 Proof of Theorem 13.1 



The proof is divided into four steps 

4.1 Galerkin approximation 

Let us define 

H m = span {exp {ik ■ x} : \k\ < m} . 

We look at the finite-dimensional equation obtained by keeping only the first m Fourrier 
modes. In order to use classical tools for systems of ordinary differential equations we need 
that / belongs to C([0,T], V' 1 ). To do so, we extend / outside [0, T] by zero and we set 

It f 
f e = Pe*f where p e {t) = -p(-), < p(s) < l,p(s) = for \s\ > 1, and / p = 1. So the 

e e Jr 
approximate sequence f e is very smooth with respect to time for p smooth and converges 

to / in the sense that 

(4.1) f e ->f strongly in L 2 ([0, T]; V^ 1 ) when e -> 0. 

The Galerkin approximation of eqs. (|2.12|) with 6 = | is given by 



(4.2) 



dii m 

— - + vP m Au m + P m B{P m vT, P m u m ) = P m fi 

Ut rn 

u m = a^Aiu" 1 + u m , 
u m (0) = P m u . 



Where for some m € N and all — 1 < s < 2, P m (w) = Yl\k\<m ™ k ex P ' x } : ^ s —> Hm 
is the orthogonal projector onto the first m Fourier modes that verifies (see in [32] for 
more details): 

(4.3) ||-fm||£(v s ,V s ) — 1; an d for all v € V s : P m v^v strongly in V s when m — > oo. 

The classical theory of ordinary differential equations implies that eqs. (|4.2I) have a unique 
C 1 solution u m for a given time interval that, a priori, depends on m, such that u m = P m u 
and V • u m = 0. Our goal is to show that the solution remains finite for all positive times, 
wich implies that T m = oo. 

In the next subsection, we find uniform energy estimates for this solution with respect to 
m. 
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4.2 Energy estimates in H 

We follow here a similar method to the one used for the Navier-Stokes equations (see |39j). 



Lemma 4.1 Let f € L 2 ([0, T], V" 1 ) and u € H, there exists K X {T) and K 2 (T) inde- 
pendent of m such that the solution u m to the Galerkin truncation \4-%\ satisfies 
(4.4) 

II u "1l 2 ([o,t],v) ^ K i( T )' f° r allT> 0, where Kx(T) = i M|wo||lr + ^II/IIz2([ 0) t],v- 1 )) 



and 

(4.5) \\u m \\l~ i[0 , T ],H) < R 2{T), for all T > 0, where K 2 (T) = 1/(^1 (T)). 

Proof Taking the L 2 -inner product of the first equation of (|4.2j) with and integrating 
by parts, using (|2.6p . the incompressibility of the velocity field and the duality relation we 
obtain 

" .mi|2 1 ,.i|v7..mi|2 _ f td -e . ~.m j m ^ 11 jjii ii„.mi 



(4.6) || w ™||^ + z ,||Vu m ||J/ = / P m /i.u TO da;<||/|| v -i||u m || v i. 

z at J T3 m 

Using Young inequality we get 

(4-7) |ll« m ||^ + H|V W m ||^ < l -\\ff v -,. 



In particular, this leads to the estimate 

(4-8) sup \\u m \\ 2 H <\\u4 2 H + -II/" 2 

te[o,T m ) 



^IIJ" Hl2([o,t],v- 1 ) 



This implies that T m = T. Indeed, consider [0,T™ ax ) the maximal interval of existence. 
Either T™ ax = T and we are done, or T™ ax < T and we have lim sup t _>. r Tma x\- \\u m (t)\\ 2 H = 
00 a contradiction to (|4.8j) . Hence we have global existence of u m for / € L 2 ([0, 00), V~ x ) 
and Uo £ H and hereafter we take an arbitrary interval [0, T] and we assume that 
feL 2 ([0,T],V~ l ). 

Integrating (|4,7p with respect to time gives the desired estimate (|4.4p for all t € [0,7"]. 
Now we use the regularization effect described in Lemma f2.2l to get the following estimates 



Lemma 4.2 Under the notation of Lemma\4- 1\ we have 



(4.9) ||w m || 2 3 <(-) Ki(T), for all T > 0, 



anc 



(4.10) llt^lr ! < - LfiT 3 (r), /or o/Z T > 0. 

By interpolation (see Lemma [24]) between L 2 ([0,T],V2) and L°°([0, T], V?), we can 
deduce the following result 
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Corollary 4.1 Under the notation of Lemma \4-1\ we deduce 

(4.11) \\u m \\l H[0 ,T],v) < (^) KxpfftK^T) 1 / 2 , for all T > 0. 

Our next result provides an estimate on the time derivative of u rn . 

Lemma 4.3 Let f G L 2 ([0, T], V~ l ) and uq G H, there exists K%(T) independent of m 
such that the time derivative of the solution u rn to the Galerkin truncation {4-<ty satisfies 

du m 

( 4 - 12 ) ll^lliwLV-) ^3(T), 

where 

K 3 (T) =4(v 2 + C (1/a) K 2 {T)) K^T) + 2||/||^ 2([0iT])V _ 1) . 

Proof. Taking the ii" -1 norm of (|4.2f) . one obtains : 
dii m 

(4-13) H^T 11 ^- 1 - "IKIIvi + + ll/llv-i 

where we note that 

\\B(u m ,u m )\\ v -i = sup {(B(u m lU m ),w);w G V\\\w\\ v i < l}. 

It remains to show that ||S(w m , u m ) \\y-i is bounded. Indeed, we have from Lemma [27 
that 



(4.14) \\B{u m , tx m )|| v -i < CH^IL i\\u m \\ v . 

It follows that 

(4-15) H^llv-< 4(^ + 6-11^11^)11^^ + 211/11^. 

Integrating with respect to time, and recalling that U™ G L°°([0, T], V? ) and w m G 
L 2 ([0, T], V), it follows from Lemma 14.11 and Lemma 14.21 that 

ft du m 

(4-16) ^ II^T'lv-i <4(^ + C(l/a)^ 2 (T))A' 1 (r) + 

4.3 Passing to the limit m — > +00 

From Lemmas 14.11 and 14.31 and from the reflexivity of the appearing Banach spaces we 

du m 

can extract a subsequence of u m and — : — such that u m converge weakly to some u in 

dt 

du r ^ (I iii 

^([O,^,^ 1 ) and - — — converge weakly to some — — in L 2 ([0,T], V- 1 ) respectively. 

Now the interpolation Lemma of Lions and Magenes (Lemma I2.6H implies that u € 
C([0,T];H). 

In order to show the convergence of u m to u in C([0, T], H)PiL 2 ([0, T], V), we need to show 
that for wo £ H, the sequence u m is a Cauchy sequence in C([0,T], H) n L 2 ([0,T],V). 
We know that u m G C([0, T], if) n L 2 ([0, T], V). The difference u m+n - u m satifies 

j/ m+n _ m\ 

(417) ^ ^ '- + P m B{{u" 1+n - tf"),ti m ) + P m B(u m+n , {u m+n - u m )) 

+vA(u m+n - u m ) = 0. 
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By taking u m+n — u m as test function in (|4.17|) we get 

l^.\\ u ™+n _ u m\\2 H + v \\ u ™+n _ 

(4.18) < \{B({u m+n - u m ), u m ), u m+n - u m )\ 

+ (Pm+nf^_ ~ Pm f± , U™ +n - U™ 



v-^v 1 

Lemma 12.41 combined with Young inequality give 

\(B((u m+n - v™),^), u m+n -u m )\< -\\u m+n - IT™ III i \\u m \\ 2 v 

+^\\u m+n -u m \\ 2 v 

The duality norm combined with Young inequality give 



(4.20) 



Pm+nf — Pmf ±,U m+n — U m \ <— \\P m+n f 1 — P m f± \\ \- 

m+n m / T/" 1 \f V m+n m 



+-||w m+n - u m f v 



Thus we get 



^||u m+n - u m f H + z^||u m+n - u m \\ 2 v < Ca^u- l \\u m+n - u m f H \\u m f v 



(4.21) dt 



By Gronwall inequality we get 

\\u m+n -U m \\ 2 H 



+ — H-Pm+n .f 1 — -Pm/ J_ ||^-1 



(4.22) * (ll- m+n (°) " + HW^-^/ifv-) 



T 

xexp/ a" 1 ^" 1 !!^" 1 !!^^ 
We know that u m G L 2 ([0, T], V), thus there exists C(a, v) > such that 



exp [ a^u^W^fdt < C(a,v). 
Jo 



We observe that 

r-T 



m+n m * J Q m+n m+n ^ 



Pm + n f 1 Pmf_L 1 — / 1 1 Prn.-\-n. f 1 -P m jf 1 Hi/- — 1 



(4.23) 



T 

+ / H-Pm-f 1 — -Pmf 1 || v -1 

q m+n m 



For / we have 

||Pm+Ti.,f 1 — Pm f 1 1 



because f i G L 2 ([0, T], V^ 1 ) and 



|-Pm.+T7.f 1 — Pm f 1 < ||P m -)- n / 1 < G L 1 [0,T]. 

m+n m+n m+n v 
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Using the dominate convergence theorem, we conclude that / tends to zero when m tends 
to oo. 

Similarly by using the dominate convergence theorem combined with the fact that 

\\Pmf 1 — Pmf ± \\v-l — > 
m-\-n m 

and 

\\P m f^-P m f±\\ 2 v -i < \\f^-f±f v -i < WfWl-i e^[o,r], 

m+n m m-\-n m 

we obtain that II tends to zero when m tends to oo. 

Since uq £ H then ||w m+ ™(0) — ii m (0)||| J - converge to zero when m goes to oo. By in- 
tegrating (OTj) . we deduce that u m+n - u m tend to zero in C([0, T],H) n L 2 ([0, T], V). 
This implies that w m is a Cauchy sequence in C([0, T], i3") n L 2 ([0, T], V). 

Concerning the initial data, we can check that w(0) = Uo. Thanks to the result above we 
have u m converge to u in C([0, T]; H). in particular it m (0) converge to it(0) in i3". In the 
other hand, we have that it m (0) = P m u o an d u m (0) converge to uq in H . The unicity of 
the limit in H allows us to deduce the result. 

It is obvious from Lemma 12.21 that u m converges to u in C([0, T]; V 1 ^ 2 ). 
Moreover, 

[ {Au m ,<t>)dt = - [ (A l / 2 u m ,A 1 / 2 cj))dt — I (A 1 / 2 u,A 1 / 2 (f ) )dt= [ (Au,<j>)dt, 
Jo Jo Jo Jo 

for all </> G ^([O^.V 1 ). 

Thus, Au m converges weakly to Au in L 2 ([0, T], V^ 1 ) as m — s> oo. 

We finish by showing that the non linear term B(u m ,u m ) converges weakly to B(u,u) 
in L 2 ([0, T], V" 1 ) as m — >• oo. From the properties of the trilinear form we have 

T f-T 



f {B{u m ,u m )^)dt- [ (B(u,u),(f))dt 
o Jo ? 

<l {(Btu 171 -u,u m ),(/))\dt+ [ \{B(u,u m -u),4>)\dt 
Jo Jo 

and by using Lemma 12.41 combined with Holder inequality we get 



I \{B{u m -u,u m ),(f>)\dt < \\u m -u|| L o C([0iT]iV i/2 ) ||u m || L 2 ([0iT]iV i ) ||^|| i 2 ([0iT])V i ) , 

J 

and 

I \{B(u,u m -u),(j))\dt < ||u|| Loo([0jT]]V i/ 2) ||-u m -«|li2([o,ri,v 1 )ll^ll£ a ([o,Ti 1 v 1 )- 

J 

Thus B(u m ,u m ) converges weakly to B(u,u) in L 2 (([0, T], V^ 1 )) asm4 oo. This im- 
plies that P m B(u m ,u m ) converges weakly to B(u,u) in L 2 (([0, T], V" 1 )) asm^ oo. 
The convergence of P m f± to / in L 2 ([0,T], V' 1 ) is obvious because / G L 2 ([0,T], V' 1 ). 

m 

We have shown that u satisfies (12.120 viewed as a functional equality in V^ 1 . 
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4.4 Uniqueness 

The solution constructed above is unique. Next, we will show the continuous dependence 
of the solutions on the initial data and in particular the uniqueness. 

Let u and v be any two solutions of (|2.12p on the interval [0, T], with initial values uq £ H 
and Vq £ H, respectively. Let us denote by w = u — v and w = u — v. Then, we can 
write the evolution equation for w as an equality in V -1 given by 

(4.24) —w + vAw + B(w, u) + B(v,w) = 

We take the inner product of (I4.24p with w . Applying the Lemma 12.61 of Lions-Magenes 
and by the properties of the trilinear form we get 

(4.25) ^IMIir + f||Vw||jf + (B(w,u),w) = 

Using Lemma 12.41 combined with Young's inequality and Lemma 12.21 we obtain 

(4.26) ^\\ w f H + u \\x/ w f H < Iln^H^ii^u^ 
Using Gronwall's inequality we conclude that 

1 4 



w \\l°°([o,t\,h) ^ ll w oll/J ex P---^i- 



a v 



We have shown the continuous dependence of the solutions on the inital data in the 
L°°([0, T], H) norm. In particular, if w = then w = and the solutions are unique for 
all t € [0,T]. 



Remark 4.1 Since T > is arbitrary, the solution above may be uniquely extended for 
all time. 

Remark 4.2 The pressure is absent in equations $2.12\) once we can recover it by using 
the de Rham Theorem Generally, the existence of the pressure is not obvious and the 
pressure may not exist 

In our situation, as we work with periodic boundary conditions, we can proceed differently. 
We take the divergence of \l-4% this yields the following equation for the pressure 



(4.27) - Ap = V • (V • uu), 

where uu is the tensor (TPu^Kij^. 

We recall that for any divergence-free field u, we have 

(u ■ V)u = V • (uu). 

Since (u ■ V)u 6 L 2 ([0,T], H^ 1 ) . We conclude from the classical elliptic theory that p is 
bounded in L 2 ([0, T], L 2 (T 3 )). 

Lemma 4.4 (Uniqueness of the pressure) The pressure is uniquely determined by the 
velocity. 
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Proof. Let (u,p) and (u,q) be two weak solutions of (1.1) (with 9 > j) such that 
Jr 3 P = Jt 3 9 = 0- Then for every € iJ" 1 and for almost every t € M + we have 

(V(p-g),0) = O, 

combined with the condition / p = q = gives p = q. 

JT 3 ila 



This concludes the proof of Theorem 13.1 

5 Proof of Theorem ET21 

5.1 Energy estimates in V? 

Let us write the energy inequality in the space Va. Taking the scalar product of eqs. 
T2l) with u m it turns out, due to the divergence-free condition, that 



(5.1) l^ii^ni +-11^^11^ = c^^(^,^),^)^, 

By definition of the scalar product on Va, we get 

(5.2) (P m B(U™,u™),«™) j < ||J9C« ra ,« m )ll„-* l|V« ra ||„ 4 



The Sobolev embedding in Theorem 12.11 together with a Holder estimate give 

\{B{u"\u m ),\V\u m ) L 2\ < c\\B(u m ,u m )\\ v - 1/2 \\Vu m \\ vl/ 2 

/ r o\ < cllfiCu^.U^HiS/allVu^llyi/a 

v / nllsrtnll llY7i.™il IIV7/».wi|l 

S C\\U \\l 6 \\vU \\L 2 \\vU ||yl/2 
S C||lt ||V|| U l|V||VW ||yl/2. 

By the interpolation inequality between V 1 ^ 2 and V 3//2 we get 

(5.4) \(B(u m , u m ), \ V\u m ) L 2 1 < llu^Hvll^ll^/allVti^ll^. 
Using the fact that 

(5.5) (P m f±,u m )vV* ^ 11/11 v-V 2 ||Vw m || vl/2 
combined with the Young inequality, we infer 

(5.6) J^WW + ^l|V^|| 2 vl/2 < ^3 llt^llvll^llvi^ + ^ll/Hv-vs- 
The following Lemmas are crucial to the proof of Theorem 1.1. 

Lemma 5.1 Let f £ L 2 ([0, T], V~^) and u 6 V^, there exists Mi(T) and M 2 (T) inde- 
pendent 0/ m stic/i i/iai t/ie solution u m to the Galerkin truncation \4-%fy satisfies 

(5-7) W um \\lmo,T],v^) ^ M i( T )' /° r aU T ^ °> 



and 



(5-8) ll um HL~([o,nvV 2 ) - M2 ( T )' /or a " * - °- 
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Proof. We proceed by two steps: 

First step: From t)5.6[) and by using GronwalFs inequality we conclude that 



(5-9) || W -||^ ([0)T]jVl/2) < M 2 (T), 

where M 2 (T) is given by 

m 2 (t) = (kh 2 v1/2 + ^ii/iii W] ,v-v 2) ) ex p^^cn 1/2 ^( T ) 1/2 - 

Second step: Integrate the original inequality (|5.6p on [0,T]. We get 
(5l0) 

l|w m (T >a; )||^ 1/2 +z/^ ||Vt* TO (t,a;)||^ 1/a d< < \\u \\ 2 vl/2 + l\\f\\l 2{[0>nv - 1/2) 

+-J^M 2 (T)K 1 (T) l / 2 K 2 {T) 1 / 2 . 

We set 

M 1 (T) = \\u \\ 2 vl/2 + ^||/||i 2([0jT])V - 1/2) + ^M^K^Tf^K^Tfl 2 . 
Thus -u m e L°°([0,T], V 1 / 2 ) nL 2 ([0,T],F 3 / 2 ) for all T > 0. 
We deduce from Lemma 12.21 that 

|2 

Il 2 ([o,t],v 2 ) 



(5.11) ll^ m |l?2anri < (l/«)Mi(T), for all T > 



and 

(5.12) ^^-([o^^i) < (l/a)M 2 (T), for all T > 

Our next result provides an estimate on the time derivative of u m . 

Lemma 5.2 Let f € L 2 ([0,T], V~ 1/2 ) and u e V 1/2 , there exists M 3 (T) independent of 
m such that the time derivative of the solution u m to the Galerkin truncation ^ -S\ satisfies 

du m 

( 5 - 13 ) w^r\\U[o,T],v- i/2 ) - Ms(T) ' for al1 T - °' 

where 



(5.14) M 3 (T) = (4(^ + C7M 1 (T))M 2 (T) + 2||/||^ M)V _ 1/a; 

Proof. Taking the V^ 1 ^ 2 norm of (14. 2b . one obtains: 
du m 

(5.15) H-^-llv-Va < "ll«1vV» + \\B(v*»,u m )\\ v - yi + ||/|| v -i /a 
where we note that 

\\B{u"\u m )\\ v - 1/2 = sup {(Bin™, u m ),w);w^V l/ \ IMI1/2 < l}- 

It remain to show that 1^(11™, «' m )|| v -i/2 is bounded. 
Indeed, we have from Lemma 12.41 that 
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(z.-ta\ ll-^V" )« ;ilv~ 1 /2 S 1 1 ^1/2 1 1 -*X || v 3/2 

It follows that 

(5-17) W^fWl-^ <4(^ + C||^||^/ 2 )||« m || 2 vl/2 +2||/|| 2 v _ 1/2 . 

Integrating with respect to time, it follows from Lemma 15. II that 
r T du m 

(5-18) J H-^-llv-Va <4^ 2 + CM 1 (r))M 2 (T) + 2||/||2 2([oT]iV/ _ 1/2) 

yields (|5,13p and the proof is complete. 



5.2 Passing to the limit 

We are now ready to take limits along subsequences of u m to show the existence of solutions 
that satisfy (|2.12p viewed as a functional equality in V~z . 

From Lemmas 15.11 and 15.21 combined with Alaoglu compactness theorem we can extract a 

du m 3 du 

subsequence of u m and — - — that converge weakly to some u in L 2 ([0,T],V2) and — m 

at at 

L 2 ([0,T],V~^) respectively. 

Now the Lions-Magenes Lemma 12.61 implies that u £ C([0, T]; V2) and this subsequence 
u m converge to u in C([0,T]; H) and in particular we have u(0) = Uq. Thus is obvious 
that u" 1 converge to u in C([0, T]; V2). 

3 3 

Fix e such that < e < 1/2, Since is compactly embedded in V^~ <L then Theorem 



21 implies that there exists a further subsequence of u m , denoted by u m \ that converges 
strongly in L 2 ([0,T];V^~ £ ). Moreover, 

f 1 T r .C ^Z" 1 

(Au mi ,<f))dt = [ {u m \Acj))dt^ [ (u,A(f>)dt= [ (Au,<f>)dt, 
Jo Jo Jo 

for all 4> G L 2 ([0,T],V^). Thus, Au m converge weakly to Au in L 2 ([0, T], V~^) as 
i — > 00. 

We finish by showing that the non linear term B(H mi , u mi ) converge weakly to B(u,u) 
in L 2 ([0,T], V~z) as i ^ 00. From the properties of the trilinear form we have 
rT pT 

(B(u mi ,u mi ),<j>)dt- / ( 
Jo 

rT r T 



B(u" l *,u ,n *),^)cft- f (B{u,u),(/))dt 

r t 
<[ \(B(u rrii —u,u mi ),(j))\dt + [ \(B(u,u m > -u),<j))\dt 
Jo Jo 



by using Holder inequality combined with Sobolev injection Theorem and Poincare in- 
equality we get 

C \(B(u m > -u,u m *),(j))\dt+ [ T \(B(u,u mi -u),cf>)\dt 
Jo Jo 

<[ \\u m '-u\\ i\\u mi \\ 3 U\\ idt+ [ \\u\\ i + \\u m *-u\\ 3 U\\ idt 
< \\6\\ 1 \\u mi -wll 1 llii mi ll 3 

_ llYll L 2 ([0,T],V?) 11 "L°°([0,T],Viy "L 2 ([0,T],V2") 

+cJ|<£|| „ 1 \\u II rooz-m r\ tzIi lit*" 1 * — till „ 3 . . 
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Thus B(u nii ,u mi ) converge weakly to B(u,u) in L 2 ([0,T], V 1 ^ 2 ) as i — > oo. 
We have shown that u satisfies (I2.12D viewed as a functionnal equality in V^ 1 ^ 2 . 



5.3 Uniqueness 

The solution constructed above is unique. Next, we will show the continuous dependence 
of the solutions on the initial data and in particular the uniqueness. 
Let u and v any two solutions of (jl.4l) on the interval [0, T], with initial values Uq € V? 
and vq G V^, respectively. Let us denote by w = u — v and w = u — v. Then from 
(12.121) we can write the evolution equation for w as an equality in V~2 given by 

(5.19) ^j-w + vAw + B(w, u) + B(v, w) = 

We take the inner product of I5.19l with w, applying the Lemma of Lions-Magenes I2U1 and 
by the bilinearity of B we have 

(5.20) + HI Viu + (B(w,u),\V\w) L 2 + (B(v,w),\V\w) L 2 =0 
The first non linear term is estimated by 

( 521 ) \{B(w,u),\V\w) L2 \ < ||wV«|| v -i/2||V«>|| v i/2 



< C||«7|| v l/2 ||Vu|| v l/2 ||Vt«|| v l/2, 



and by Young's inequality, we obtain 



(5.22) \{B(w,u), \V\w) L 2\ < -||w||^i/ a ||« ||^3/2 + ^||Vw||^ 1/2 . 
The second non linear term is estimated by 

\(B(v,w),\V\w) L 2\ < \\vVw\\h\\Vw\\h 

(5.23) < CllwIlvllVtollyi^llVtollff 



and by Young's inequality, we obtain 



, n 3/2 n , 1 1/2 

< C\\v\\v\\w\\^^\\w\\^ !2 



(5.24) \(B(v,w), \V\w) L 2\ < -||u||^|H| 2 + ^-\\Vw 



,.yl/2 -T 4 || V U/|| v i /2 . 



From the above inequalities we get 

d 1 1 , | o I, , , 1 1 o C 



( 5 - 25 ) j t \\™\\vV*+ v \\ Vw \\W ^ -W W WW (ll»llv + ll«llv3 /3 ) 

Using Gronwall's inequality we conclude that 



w 



2 ^ II... ||2 



2 , 1/2 exp^^y + J 



L tx) ([0,T],V 1/2 ) — II Oil v 1 / 2 ^ ~jj \ / II^HV^T" / ll"llv 3 / 2 



Note that the above estimate does not depend on a, thus we have not used the regular- 
ization effect to get the uniqueness in V"2 and this is not surprising because the solutions 
of the Navier-Stokes equations with initial data in V? are unique. 

We have shown the continuous dependence of the solutions on the inital data in the 
L°°([0, T], V 1 / 2 ) norm. In particular, if w = then w = and solutions are unique for 
all t S [0, T]. Since T > is arbitrary this solution may be uniquely extended for all time. 
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Remark 5.1 The pressure is absent in eqs. once we can reconstruct it from u and 

u (up to a constant), if necessary, see Lemma {4~4\ and Remark \4-2\ We take the divergence 
of Jl.^p . This yields the following equation for the pressure 

(5.26) - Ap = V • (V • uu). 

Where uu is the tensor (u t u :i )i<ij<s. 

We recall that for any field u with divergence is equal to zero, we have 

(u • V)w = V • (uu). 

Since (u-V)u £ L 2 ([0, T], H~ l l 2 ). One concludes from the classical elliptic theory that p 
is bounded in L 2 ([0,T],H 1/2 {T 3 )). 



This finish the proof of Theorem 13.21 

6 Relation between NSE and Leray-a: construction of suit- 
able solutions 

The regularized solution contructed above with 6 > | is unique. So the solution is suitable. 
In this section we will construct a suitable weak solution to the Navier-Stokes equations 
by taking the limit when a tends to zero to the regularized solution. 
We start by the following Lemma: 

Lemma 6.1 Let (u a ,p a ) be the unique solution of |L^| ) with 6 > \ then (u a ,p a ) verifies 
the following local inequality 



T 

Jj 3 



2u I I \X7u a \ 2 <fi dxdt 



(6.1) = / / |« a | 2 (4>t + vAcf)) + (|w a | 2 w a + 2pu a ) ■ V<fi dxdt 

Jo Jt 3 
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rT 



+2 / fu4> dxdt, 
Jo Jt s 

for all T € (0, +oo] and for all non negative fonction <p € C°° and supp <p CC T3 x (0, T). 

Proof. We take 2(pu a as test function in (|1.4jl . We note that the condition 9 > 1/4 ensure 
that all the terns are well defined. In particular the integral 

2 / i^Vu a • u a (f) dxdt 
Jo Jt 3 

is finite by using the fact that u^Vu a € L 2 ([0, T];H~ l ) and 2(j)U a G L 2 ([0, T];H l ). 
An integration by part combined with (j>(T, ■) = 0(0, •) = and the following identity 



(6.2) 2 / u a \7u a ■ u a (j) dx = j u a \u a \ -V/fidx 

Jr 3 

yield that (u a ,p a ) verifies 



In order to take the limits a — > over (|6.ip . we need first to show that for all u a £ 
LP(0,T;L p (T 3 ) 3 ), 2<p< 10/3, we have: 

(6.3) u^ -> u strongly in L p (0, T; L p (T 3 f) for all 2 < p < 10/3. 

This is the aim of the two following Lemmas. 
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Lemma 6.2 Let 9 € R+, < (5 < 29, s G R and assume taai u £ H s . Then u € -H" s+/3 
and 

(6-4) INff.+^ < ^gll«llff- 



Proof. When u = exp {ifc • a;}, then 

fee/ 

(6.5) ^ = E r+j%F exp{ik-x} 

kei ' ' 

Formula (16. 5h easily yields the estimate 



|u||^ < (^P (1 + ^l| fc | 2e)2 ) INI"* 

2^ 



(6.6) . / 1 ( «W_VVll lir 



a 2/3 fee/ Vi + a^W 29 

1 II lir 

where we have used the fact that 

(6 - 7) a d+^i*p J sl>a " ,SJS * 

Lemma 6.3 Assume u a belongs to the enegy space of solutions of the Navier-Stokes 
equations, then 

o 10 

(6.8) iZ^ -> it strongly in L p (0, T; L p (T 3 f) for all 2 < p < y . 

3p-6 

Proof. From the Sobolev injection H 2 ? '—t L P (T 3 ) , it is sumcent to show that 

(6.9) / || || 3p _6 — 0, when a — > 0. 

io if^P~ 

From the relation between and u a we have 

(6.10) \\u^-u a \\ p 3p _ 6 < a 26lp ||TZ^|| p 3p _6 



Y28 



Lemma 16.21 implies that 



10-3p 

(6.11) / \\u^- u a \\ p 3p _ 6 dt < oT^r 2 - I \\u a \\ p z dt. 

JO H~^p~ JO Hp 

Recall that 

(6.12) / ||w Q || p 2 dt < oo fo any 2 < p < oo. 

Jo hp 

This yields the desired result for any 2 < p < 10/3. 

In order to show that (u a ,p a ) gives rise to a suitable solution of the Navier-Stokes equa- 
tions, it is necessary to take the limit a — > 0. We have the following theorem: 
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Theorem 6.1 Let (u a ,p a ) be the solution of Ji,^| j ; in the a — s> limit, (u a ,p a ) tends to 
a suitable solution to the Navier-Stokes equations. 



Proof: By a classical compactness argument [13] . we deduce that u a approaches u 
strongly in L 2 ([0,T], H) for all T > when a tends to zero, where u is a weak solu- 
tion to the Navier-Stokes equations . It remains to show that (u,p) verifies the following 
local energy inequality 



2u I [ \Vu\ 2 4> dxdt < [ [ \u\ 2 (4> t + vA(f>) + (\u\ 2 u + 2pu) • V</> dxdt 
(6.13) ^° ^ T3 ^° ^ T3 



1 3 J13 

fT 



+2 / / fu<p dxdt. 
'o Jt 3 



To do so, we need to find estimates that are independent from a. Using the fact that 
(u a ,p a ) belong to the energy space: L°°([0, T]; H) n L 2 ([0, T]; V) n if ([0, T}; if (T 3 ) 3 ) 
and from the Aubin-Lions compactness Lemma (the same arguments as in section 4) we 



can find a not relabeled subsequence (u a ,p a ) and (u,p) such that when a tends to zero 
we have: 


(6.14) 


u a - 


x u 


weakly* in L°°([0, T];H), 


(6.15) 


u a - 


x u 


weakly in L 2 ([0, T}; V) n Lt ([0, T]; L% (T 3 ) 3 ), 


(6.16) 


Pa ~ 


±p 


weakly in L§ ([0, T]; L§ (T 3 )), 


(6.17) 


Ua ~ 


> u 


strongly in L p ([0, T]; L P (T 3 ) 3 ) for all 2 < p < — , 

3 


(6.18) 


Ua~ ~ 


> u 


strongly in L p ([0, T]; L P (T 3 ) 3 ) for all 2 < p < — , 

3 


(6.19) 


Pa ~ 


X p 


strongly in L p ([0, T]; L P (T 3 )) for all 1< p < -. 

(J 


These convergence results allow us to take the limit in ([6.1]) in order to obtain, by using 
the weak lower semicontinuty of the norm in L 2 ([0, T]; V) 






lim inf 

a^O 


/■T /-T 

/ ||Vu a ||H dt > / ||Vu||h di, 


that (it, 


p) verifies the local enerev inequalitv (|6.13p. 



7 The deconvolution case 

7.1 The modified deconvolution operator 

In this section, we introduce the modified deconvolution operator which interpolate the 
usual deconvolution operator introduced in [23] . 

Let a > 0, s > -1, < 9 < 1, u G H s and let u G H s+2e be the unique solution to the 
equations 

(7.1) a 28 (-A) e u + u = u. 

(7.2) V-u = V-TZ = 
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We also shall denote by G the operator 

^ H s + 20 — > //\ 



(7.3) 



w — >■ a 2e (-A) y i<; + («. 



Therefore, one has 

(7.4) u = G~ 1 u. 

Let us consider the operators 



Dn = £(/ 



TV 

i — 1 \ n 



™ k e ikx . 



71=0 

and 

(7.5) = D N (u). 
A straightforward calculation yields 

(7.6) ^^^.^^-(^j^)^ 

One can prove the following (see in |23^ I27j with 9 = 1): 

• Assume u G i/ s . Then ffjv(«) G ii" s+2e and ||.H'jv(u)|| jr .+» < C(iV, a)||u||fr», 
where C(N,a) blows up when a goes to zero and/or N goes to infinity. This is due 
to the fact 

/ / a 2e \k\ 2e \ N+1 \ N + l 

[ \l + a 2e \k\ 2 °) )~a™\k\™ ™ 1 |o ° 00 ' 

• The operator Hn maps continuously H s into H s and H-ff/vll/^/p^s) = 1- 

• Assume u G H s . Then the sequence (u a ) a> o converges strongly to u in the space 
H\ 

• Let u G -H" s . Then the sequence (i7/v(*0);veN converges strongly to n in when 
./V goes to infinity. 

7.2 The deconvolution model: well-posedness results 

We consider here a family of equations interpolating between the Navier-Stokes equations 
|25j and the Leray-deconvolution model [23] with periodic boundary conditions. 

— + H N (u)-Vu-uAu + Vp = f in R + x T 3 , 

V • u = V • H N (u) = 0, u= H N (u) = 0, 

h A Jt 3 

u(t,x + Lej) = u(t, x), 



(7.7) 



Where H^(u) is an interpolating deconvolution operater introduced above. 

When N = 0,we obtain the equations (|l,4p . Similarly, we obtain the same results of 
well-posedness for the interpolating model (|7.7p . 
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Theorem 7.1 For any T > 0, let f G L 2 ([0,T],V~ 1 ) and u G H. Assume that 9 = 1/4. 

Then there exists a unique solution (u,p) := (Wo!,JV)Pa,iv) ^° egs. J 7. 7| j t/iat satisfies 

nil 

u G C([0,T];2f) nL 2 ([0,T];y) and — € L 2 ([0, T]; F" 1 ) andpS L 2 ([0, T], L 2 (T 3 )). 
SW/i taat it verifies 

(^- + vAu + B(H N (u),u)- f,(t>) = 

\ at / v-^v 1 

for every (p £~V and almost every t G (0, T). Moreover, this solution depends continuously 
on the initial data uq in the L? norm. 

And 

Theorem 7.2 ForanyT>0, let f G L 2 ([0,T],V~^) and u G . Assume that 9 = \. 

Then there exists a unique solution (u,p) := (u ai N,p a ,N) t° e Q s - ( f7 7p that satisfies 

i ^ fin i i 

u eC([0,T];V2)nL 2 ([0,T];V2), — G L 2 ([0,T],V-2) and p G L 2 ([0, T], H 5 (T 3 )). 

SW/i t/iat it satisfies 

/^ + vAu + B(H N (u),u)-f,4>) i =0 

\ ai / V~^,V"^ 

for every <j> G and almost every t G (0, T). Moreover, this solution depends continu- 
ously on the initial data Uq in the Hz norm. 

Since the proof of Theorem 17.11 and 17.21 are similar to that of Theorem 13.11 and 13.21 we omit 
their proof here. 



7.3 Convergence to a suitable weak solution to the NSE 

By a similar argument to that in [23], it is possible to prove that up to a subsequence 
the solution for the deconvolution model (with = -i) converges when a — > or/and 
N — > +co to a weak Leray solution to the Navier-Stokes equations. Next, we will show 
that the deconvolution regularization gives rise to a suitable weak solution to the Navier- 
Stokes equations. We need first 

Lemma 7.1 Let {u aj N,Pa,N) be the unique solution of J7. 7| ) with 9 > | then (u a! N,Pa,N) 
verifies the following local inequality 



(7.8) 



2v I I \X7u a ^\ 2 4> dxdt = / / \u a ^\ 2 (<j) t + i>A<j)) dxdt 
Jo Jt 3 Jo Jt 3 



+ / / (\u a , N \ 2 H N (u a . N ) + 2p a>N u a , N ) ■ V4> dxdt + 2 / / fu a , N (f> dxdt, 
Jo JTs Jo Jt 3 



for all T G (0, +oo] and for all non negative fonction (f> G C°° and supp (f> CC T3 X (0, T). 
Proof. See Lemma loTTl 

Then we have 
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Lemma 7.2 Assume u a ^ belong to the enegy space of the solutions of the Navier-Stokes 
equations, then 

, 10 

(7.9) H N (u a , N ) -> u strongly in L p (0, T; L P (T 3 ) ) for all 2 < p < — . 

3 

3p— 6 „ 

Proof. From the sobolev injection H 2 p ^-?> L P (T 3 ) , it is sufficent to show that 

(7.10) / \\HM{u as N) — u a ^\\ p 3p _ 6 dt — > 0, when a — > or/and iV — >• oo. 
Jo h^p~ 

From the relation between H^(u a) N) an d u Q) jv we have 



a^e/ aV 

(7.11) fceJ 



/ 20ii 1 2(9 \ 2(JV+1)\ \ 2 

, 3p-io / or* fen x x x 



- (^ ( |fc| ' [tt^w) ) ) ||u »;» 

Using the fact that l m c — >■ Z°° for any m>lwe get for p < ^ that 
(7.12) 

ii ff r ^ hp < /V ( 1 f ^ 2m(7V+lA X 1 



Now, we choose m such that m 10 p 3p > 3, and recall that 

r 

(7.13) / \\u a n\\ p 2 dt < oo fo any 2 < p < oo, 

Jo ' hp 

and 

(7 - 14) i+a^ifcr or/and S [tt^w J - L 

Using the dominate convergence theorem, we conclude the desired result for any 2 < p < 
10/3. 



J J U HP 



Remark 7.1 Note that the above result holds true forp = 10/3. 

The above L p convergence combined with the fact that u a ^ belong to the energy space of 
solutions of the Navier-Stokes equations and the Aubin-Lions compactness Lemma allow 
us to take the limit a —¥ or/and N — y oo in (|7.8p and to deduce the following Theorem: 

Theorem 7.3 Let {u a ^,Pa,N) be the solution of J 7. 7| ), Then when a — > and/or N — > 
o°; (u a ,N,Pa,N) tends to a suitable weak solution (u,p) to the Navier-Stokes equations. 
The convergence to u is weak in L 2 ([0,T];^) n L~ ([0, T]; L~ (T3) 3 ) and strong in 

5 3 

■ 



L q ([0,T];L q (T 3 ) 3 ) for all2 < q < ^ and the convergence topis strong in L q ([0, T];L q (T 3 )) 
for all 1 < q < f and weak in Ls([0,T];L3 (T 3 )). 
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8 The MHD Case 



In this section, we consider a deconvolution-type regularization of the magneto-hydrodynamic 
(MHD) equations, given by 

(8.1a) d t u - utAu + (H N (u) ■ V)u - (H N (B) • V)S + Vp + ^V\B\ 2 = 0, 

(8.1b) d t B - v 2 AB + (H N (u) ■ V)B - {H N {B) ■ V)u = 0, 

(8.1c) V • B = V ■ H N (B) = V • u = V • H N {u) = 0, 

(8.1d) B{0) = B , u(0) = u , 

where the boundary conditions are taken to be periodic, and we also assume as before 



that / udx = Bdx = 0. 

Jt 3 Jt 3 

Here, the unknowns are the fluid velocity field u(t,x), the fluid pressure p(t,x), and the 
magnetic field B(t,x). Note that when a = 0, we formally retrieve the MHD equations 
and the MHD-Leray-a equations are obtained when N = and 9 = 1. 
Existence and uniqueness results for MHD equations are established by G. Duvaut and 
J.L. Lions in |llj . These results are completed by M. Sermange and R. Temam in [31] . 
They showed that the classical properities of the Navier-Stokes equations can be extended 
to the MHD system. The use of Leray-a regularization to the MHD equations has received 
many studies see [28]. The idea to use the deconvolution operator from [23] to MHD equa- 
tions is a new feature for the present work. 



8.1 Existence, unicity and convergence results 

First, we establish the global existence and uniqueness of solutions for the MHD-Deconvolution 
equations ([%7T]) for 9 = 1/4. 
We have the following theorem: 

Theorem 8.1 For 9 = 1/4. Assume u G H and Bq € H. Then for any T > 0, \8.1\) has 

a unique regular solution (u,B,p) := (u a ^ ,B aj N ,p a ,N) such that, u,B € C((0, T],H) fl 
L 2 ([0, T], V 1 ), and p € L 2 ([0, T], L 2 (T 3 )).' 
Furthermore the solution verifies 

/ [ (\Vu\ 2 + \VB\ 2 )(f> dxdt 

!> 

T 



f f \u\ 2 (fo + ^Acf)) + \B\ 2 ((j) t + u 2 A(j)) dxdt 
Jo Jt 3 

rT r 



(8.2) ■><> ;/a 

+ / / ((\u\ 2 + \B\ 2 ) H N (u) + 2pu) -V0 dxdt 

JO Jt- a 

-4 f [ (uB) H N (B) ■ V0 dxdt 
Jo Jt 3 

for all T € (0, +oo] and for all non negative fonction <f> € C°° and supp <p CC T3 x (0, T). 

Proof. We only sketch the proof since is similar to the Navier-Stokes equations case. The 
proof is obtained by taking the inner product of ()8.1ap with u, (|8.1b|) with B and then 
adding them, the existence of a solution to Problem (|8.ip can be derived thanks to the 
Galerkin method. Notice that u,B satisfy the following estimates 

(8.3) ~ {\\u(t,x)\\ 2 H + \\B(t,x)\\ 2 H ) +mm(u 1 ,u 2 ) (HVti^aj)!^ + ||VB(t, x)\\ 2 H ) < 
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The pressure p is reconstructed from it, Hn(u),B, Hjy(B) (as we work with periodic bound- 
ary conditions) and its regularity results from the fact that (Hjy(u) ■ V)w, (Hn(B) ■ V)£> € 
L^TIH- 1 ). 

It remains to prove the uniqueness. Let (ui,Bi,p±) and (u2,B 2 ,P2), be two solutions, 
5u = u 2 — u\,6B = B 2 — Bi, 5p = P2 — Pi- Then one has 

(8.4a) 

d t 8u + (H N {ui)V)5u - (H N {Bi)V)5B - uxAdu + V5p = -{H N (5u)V)u 2 + (H N (5B)V)B 2 , 
(8.4b) 

d t SB + {H N (ux)V)6B - (H N (Bx)V)8u - u 2 A5B = (H N (5B)V)u 2 - (H N (5u)V)B 2 , 

and Su = 0, 5B = at initial time. One can take 5u € L°°((0,T],H) nL 2 ([0,T],V) as 
test in (j53aD and 6B € L°°((0, T],H) n L 2 ([0,T], V) as test in lpb| . Since H N (ui) is 
divergence- free, one has 



5.5a) f [ {H N (u 1 )V)5u.5u = 0, 

Jo Jr 3 

i.5b) [ [ (H N ( Ul )V)5B.6B = 0. 

Jo Jt 3 



Therefore, 



(8.6) 



and 



I \5u\ 2 + v 1 j \V5u\ 2 -[ {H N {B{)V)5B.5u 

J"., Jr* 

(H N (5u)V)u 2 .5u+ [ (H N (5B)V)B 2 .5u, 
Jt 3 



^- [ \5B\ 2 + v 2 [ \V5B\ 2 - [ {H N {B 1 )V)6u.6B 

(8.7) Mt Jf z Jt z Jt z 

= / {H N {SB)V)u 2 .5B - / (H N (5u)V)B 2 .SB. 
Jt 3 Jt 3 

One has by a integration by parts, 

(8.8) - I (H N {Bx)V)8B.5u = [ {H N (Bi)V)5u.8B. 

Jr 3 

One has by adding ([876]) . (f87T|) and using ([878]) 



_d_ 

2dt 



[ \5u\ 2 + ^-f \8B\ 2 + Vl f \V5u\ 2 + v 2 [ \V5B\ 2 
Ji 3 *<*t Jt 3 Jt 3 Jt 3 

(8.9) =- / (H N (Su)V)u 2 .5u+ [ (H N (5B)V)B 2 .5u 

Jr 3 Jt 3 

+ / (H N (SB)V)u 2 .5B - / (H N (5u)V)B 2 .8B. 

Jt 3 Jt 3 
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One has by a integration by parts, 



inequality combined with Sobolev injection 



(8.10a) 
(8.10b) 
(8.10c) 
(8.10d) 
By Young's inequality, 



(8.11a) 
(8.11b) 
(8.11c) 
(8.11d) 
By Holder 

(8.12a) 

(8.12b) 

(8.12c) 

(8.12d) 
Hence, 

(8.13) 

Therefore, 



{H^(5u)V)u 2 .5u 



(H N (SB)V)B 2 .5u 



(H N (5B)V)u 2 .5B 



(H N (5u)V)B 2 .5B 



Hn(5u) (g> u 2 : VSu, 
H N {5B)®B 2 : V5u, 
H N {5B) ® u 2 : V5B, 
H N {5u)®B 2 : V5B. 



f (H N (5u)V)u 2 .5u\<^ [ \V5u\ 2 + — [ \H N (5u)\ 2 \u 2 \ 

T 3 4 JT 3 y l JT 3 

f (H N (5B)V)B 2 .5u\<^ [ \s/5u\ 2 + — [ \H N (5B)\ 2 \B 2 \ 

T 3 4 JT 3 v l JT 3 

[ (H N (5B)V)u 2 .5B\<^ [ \V5B\ 2 + —[ \H N {5B)\ 2 \u 2 \ 

Jf 3 4 JT 3 v 2 H 3 

[ (H N (5u)V)B 2 .5B\<^ [ \VSB\ 2 + -[ \H N (5u)\ 2 \B 2 \ 
Jt 3 4 Jt 3 u 2 Jt 3 



1 

Vl J J. 



1 



1 

v 2 Jj 



\H N {5u)\ 2 \u 2 \ 2 < 



< 



\H N (5B)\ 2 \B 2 \ 



\H N {5B)\ 2 \u 2 \ 



- I \H N (5u)\ 2 \B 

v 2 Jj 3 



1 

V\ 

1 

" V\ 
1 

< — 
1 

< — 
1 

< — 
1 

< — 

^2 
1 

< 
1 

< — 

V2 



2 ^ — 



H N (5u)\\ 2 L3 \\u 2 \\ 2 L6 



H N (5u)\\ 2 i\\u 2 \\ 2 v , 



2 



H N {5B)\\i 3 \\B 2 
H N (6B)\\ 2 AB 2 \\ 2 V , 
H N (5B)\\ 2 L3 \\u 2 \\ 2 L6 
H N (5B)\\ 2 vh \\u 2 \\ 2 v , 
H N (5u)\\ 2 L3 \\B 2 \\ 2 L6 
H N (5u)\\ 2 AB 2 \\ 2 V . 



min K* 2 ) \^ 5u Kh + W H ^ B K\) («V + ^ 
< m f n ( ^ 2) (W^fv + m?v) {\\5uf H + \\SBf H ) . 



A (ii^i^ + ii^n^ < c{t) + \\SB\\ 2 H ) , 
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where C(t) = V ' , (\\u 2 \\v + llftllv) € l1 ([0,T]). We conclude that 6u = 5B = 

nun (1/1,1/2) 

thanks to Gronwall's Lemma. 

In order to deduce the local energy equality (|8.2p . we multiply the equation (|8.1a|) with 
2ucj), the equation fjS.lbf) with 2B(j), for all T € (0, +00] and for all non negative fonction 
4> € C°° with supp <j) CC T3 x (0,T), and then adding them. The rest can be done in 
exactly way as in Lemma l6.1[ so we omit the details. 

With smooth initial data we may also prove the following theorem 

Theorem 8.2 For 6 = 1/4. Assume u £ V? and Bo € V*. Then for any T > 
0, ( OOP has a unique regular solution (u,B,p) := (u at N,Ba,N,Pa,N), such that u,B £ 
C([0,T],yi)nL 2 ([0,T],vf), andp€L\[0,T\,H*(Ts))- 

Proof. We only sketch the proof since is similar to the one for Navier-Stokes equations. 
The proof is obtained by taking the inner product of (IHTTaD with |V|u, (l8~Tb) with \V\B 
and then adding them, the existence of a solution to Problem (|8.1|) can be derived thanks 
to the Galerkin method. Notice that u,B satisfy the following estimates 

H (\\ U W 2 1 + II^H 2 +min(p u u 2 ) (\\Vu\\ 2 ! + \\VBf v 12 ) 
2d,t\ YJ "vhJ \ T5 v J 

[ H N (B)VB\V\udx 

■/T 3 

f H N {B)Vu\V\Bdx . 
Jt 3 

The first non linear term is estimated by 



.14) 



< 



H^s{u)Wu\V\udx 



+ 



+ 



H N (u)VB a > N \V\B a ' N dx 



+ 



J 

Jj 3 



H^{u)Wu\V\udx 



5.15) 



< 

< 
< 
< 

< 



H N (u)Vu\\ v _i ||Vw|| v 

Hn(u)Vu\\ 3 IIVwII 1 
Hn(u) \\l6 II Vw||i,2 1| Vu 
H N {u)\\ v \\Vu 

H N (u)\\ v \\u 



1 

jf||Vtt|| 1 

-nil llyj 

1 3 

2 1 l|V-"|| 2 !. 
V? 



Where we haved used the Sobolev embedding in Theorem l2.1l together with Holder estimate 

1 3 

and the interpolation inequality between V2 and V2. 
Similarly, we can estimat the third term by 



1.16) 



H N (u)VB\V\Bdx 



< 



H N (u) 



\BP 1 ||Vtff J. 



The second non-linear term is estimated by 



H N (B)VB\V\udx 



.17) 



< 

< 
< 
< 
< 



\H N (B)VB\ 



v 



1 Vu 

2 11 

IVu 

II 6 '" 

\\H n (B)\\v\\VB\\ h 
C(a,N)\\BL,i\\B 



1 



|^(S)V5|| l 3 

\H N (B)\\ L e\\VB\\ L 2 



1 



1 

1 V? 

\Vu\\ 1 
Vu|| 1 
"Vw|| !. 

V3 
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By the same way we obtain 



H N (B)Vu\V\Bdx 



< C(a,N)\\B\\ vi \\u\\ v \\Vu\\ vi . 



(8.18) 
Therefore, 

( 8 - 19 ) <||^(«)||v||u||L||V«||L + \\H N (u)\\ v \\B\\K\\VB\\K 

+C7(a,JV)||B|| v i||B|| v ||Vu||^ 4 + C(a, N)\\B\\ vi ||tt|| v l|Vu||^ 

By Young inequality, 
d 



V \<h + llSll^j +min(, 1 ,, 2 ) g|V«ir vi + ||VB||^ 

(8.20) <— 7^ — xPM«)ll 4 HI; U + — x — dl^Nll^llsf i 

^ ^ mm (1/1,1/2) ^2 mm (1/1,1/2) ^ 2 



+ Ci ^ N \ \\B\r,\\Bf v+ C( °> N \ \\B\\* 
min(i/i,i/ 2 ) V2 11 v mm(z/i,z/ 2 ) V2 



Therefore, 



^ (|M| 2 ! + ||£|| 2 +min(v 1 ,v 2 ) (\\Vu\\ 2 , + ||V£|| 2 x 

<||«||2 f C II ff r«M|4 1 gW^Q iiBiia 1 C(a,N) 

S O 1 — : — 7 iFM tt V"i : — 7 T " V "i : — 7 \ u 

(8.21) v ^ \nim(vi,v 2 ) m.m(u 1 ,u 2 ) mm (1/1,1/2) 



h — : — 7- — -HFivHlliHMI 2 1 

mm (1/1,1*2) V2 
< C(t) fll^H 2 1 + \\B\\ 2 



||».|2 , „..„2 \ r l, 



where C(t) = — -JH N (u)\\tr + .V y J |g||^ + .V M u|lfr E ^([0,T]). 

\ mm (1/1,1/2) mm (1/1,1/2) min (1/1,1/2) / 

We conclude that u, B E L°°((0, T], V^) n L 2 ([0, T], vl) thanks to Gronwall's Lemma. 

The pressure p is reconstructed from u, Hn(u),B, H^(B) ( as we are working with periodic 
boundary conditions) and its regularity results from the fact that (Hn(u) ■ V)u, (Hjy(B) • 
V)BeL 2 ([0,T],H~ l 2). 

It remains to prove the uniqueness. Let (ui,B±,pi) and (u 2 , B 2 ,p 2 ) be two solutions, 
8u = u 2 — ui,5B = B 2 — B\, Sp = p 2 — p\. Then one has 

d t 5u + (H N ( Ul )V)Su - (H N (Bi)V)SB - uASu + V5p 
[ ' = -{H N {5u)V)u 2 + (H N (5B)V)B 2 , 



d t 5B + (H N ( Ul )V)SB - (H N (B!)V)6u - uASB 
1 ' = (H N (5B)V)u 2 - {H N {5u)V)B 2 , 

and<5u = 0, 5B = at initial time. One can take \V\Su £ L°°((0,T], V~^)nL 2 ([0,T], W) 
as test in (18^21) and |V|<fiB € L°°((0, T], V~5) n L 2 ([0, T], Vs) as test in (^231) . 
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Once we obtain by a similar way as in Theorem 18.11 and Theorem 13.21 that 

(||H| 2 i + ll^ll 2 < C(t) (\\5u\\ 2 , + \\5B\\ 2 , 
2dt V V? TV V "V2 V2 

where C(i) S L 1 ([0,T]). We conclude that Su = 5B = thanks to Gronwall's Lemma. 

Next, we will deduce that the deconvolution regularization give rise to a suitable weak 
solution to the MHD equations. 

Theorem 8.3 Let {u a ^, &a,N,Pa,N) be the solution of |7. 7\ ), Then when a — > and/or 
N — > oo, (u a n, B a N, ,Pa,N) tends to a weak solution (u, B,p) to the MHD equations. The 
convergence to u and the convergence to B are weak in L 2 ([0, T]; V) flLT ([0, T]; L~3~ (T3) 3 ) 
and strong in L q ([0, T]; L q (T^) 3 ) for all 2 < q < ^ and the convergence to p is strong in 

L q ([0,T];L q (T 3 )) for all 1 < q < | and weak in L§ ([0, T]; L§ (T 3 )). Furthermore the 
solution verifies in addition 

I [ (\Vu\ 2 + \VB\ 2 )(p dxdt 
■If 3 



.24) 



< / / |w| 2 (0i + i^i A0) + l^l 2 (c^t + i/ 2 A0) dxdt 

■JO JT 3 ^ 

+ / / {{\u\ 2 + \B\ 2 ) u + 2pu) -V(j)dxdt 
Jo Jt 3 

-4 f [ (uB) B ■ Vcf) dxdt 



JO JT 3 

for all T € (0, +00] and for all non negative fonction (f> € C°° and supp (f> CC T3 X (0, T). 
Proof. As in Lemma 17.21 we can show that 

, 10 

(8.25) H N (u a)N ) -> u strongly in L p (0, T; L p (T 3 f) for all 2 < p < — . 

O 

o 10 

(8.26) H N (B a , N ) -> B strongly in L p (0, T; L P (T 3 ) 3 ) for all 2 < p < — . 

The above L p convergence combined with the fact that u a ^N and £> Q; 7v belong to the enegy 
space of the solutions of the Navier-Stokes equations and the Aubin-Lions compactness 
Lemma allow us to take the limit a — > or/and N — > 00 in (|8.2|) . The rest can be done in 
exactly way as in [23], so we omit the details. 
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